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Spectra of the spreading layers on the neutron star surface and 
constraints on the neutron star equation of state 
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ABSTRACT 

Spectra of the spreading layers on the neutron star surface are calculated on the basis of the 
Inogamov-Sunyaev model taking into account general relativity correction to the surface grav- 
ity and considering various chemical composition of the accreting matter. Local (at a given 
latitude) spectra are similar to the X-ray burst spectra and are described by a diluted black 
body. Total spreading layer spectra are integrated accounting for the light bending, gravita- 
tional redshift, and the relativistic Doppler effect and aberration. They depend slightly on the 
inclination angle and on the luminosity. These spectra also can be fitted by a diluted black 
body with the color temperature depending mainly on a neutron star compactness. Owing to 
the fact that the flux from the spreading layer is close to the critical Eddington, we can put 
constraints on a neutron star radius without the need to know precisely the emitting region 
area or the distance to the source. The boundary layer spectra observed in the luminous low- 
mass X-ray binaries, and described by a black body of color temperature T c = 2.4 ±0.1 keV, 
restrict the neutron star radii to R = 14.8 ± 1.5 km (for a 1.4-M© star and solar composition 
of the accreting matter), which corresponds to the hard equation of state. 
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1 INTRODUCTION 

Matter accreting on to a weakly magnetized neutron star (NS) in 
low mass X-ray binaries (LMXRBs) can form an accretion disc 
which extend down to the NS surface. A boundary layer (BL) is 
formed between the accretion disc and the NS surface, where a 
rapidly rotating matter of the disc is decelerated down to the NS an- 
gular velocity. The amount of the energy, which is generated during 
this proc ess, is co mparable with the energy g enerated in the accre - 
tion disc ISunvaev & ShakurJl98ri lslbeatullin & Sunvaevll998l) . 

There is no generally accepted theory of the BL. Two dif- 
ferent approaches to the BL description are considered. First 
of them, which we will call a 'classical model', considers 
the BL between a central s tar (a white dwarf or a NS) as 
a par t of the accre tion disc jprin ^lej |l977|:jPrin^le & Savonije 



-R ~ Hbl 



(2) 



1979; Tvlenda 1981; Shakura & Sunyaev 1988; Bisnovatvi-Kogan 
1994llPopham & Naravanll995l:IPoDham & Sunvaevl200ll) . In this 
model the component of velocity normal to the accretion disc plane 
is zero. The half-thickness of the BL is determined by the same re- 
lation, as for the accretion disc: 



#bl ~ — R, 

VK 



(1) 



where c s is a sound speed in the BL and vk is the Keplerian velocity 
close to the NS surface of rad ius R. The ra dial extension of the BL 
is determined by the relation lPringlell97i 



In this classical model the accreting matter in the BL is deceler- 
ated in the accretion disc plane, along radial coordinate only, due to 
the viscosity operating within the differentially rotating BL, simi- 
larly to the accretion disc. From the observational point of view, the 
classical BL is a bright equatorial belt close to the NS surface. The 
effective temperature of the BL is higher than the maximum accre- 
tion disc effective temperature, because the BL is smaller than the 
accretion disc, while their luminosities are comparable. 

Another approach was suggested by llnoeamov & Sunvaevl 

( 1999, hereafter IS99). The BL is considered as a spreading layer 
(SL) on the NS surface. The accreting matter diffusing along the 
radial direction in the accretion disc and reaching the neutron star 
surface gains the velocity component normal to the accretion disc 
plane due to the ram pressure from the accretion disc. Then the 
matter spirals along the NS surface towards the poles. Rotating at 
the NS surface, the matter is decelerated due to a turbulent friction 
between the rapidly rotating matter and a slowly rotating NS sur- 
face. The kinetic energy of the accreting gas is mostly deposited in 
two bright belts at some latitude above and below the NS equator. 
The width and the latitude of the belts depend on the mass accre- 
tion rate. The larger the accretion rate, the wider the belts are and 
the closer they are to the NS poles. At the accretion rate close to 
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the Eddington limit (Lbl ~ i-Edd) the bright belts expand all over 
the NS surface. 

The observational difference between two BL models is not 
significant. At low accretion rates (Lbl ~ O.OlZ-Edd) the latitude 
of bright belts of the SL is small (~ few degrees) and the verti- 
cal extension of the SL is comparable to the classical BL thick- 
ness. At high accretion rates (Lbl ~ ^ Edd) the cl assical BL thick- 
ness is comparable to the NS radius (see lPopham & SunvaevfcOOll) . 
Therefore, the effective temperatures of these BL models are of the 
same order. 

In the approach by IS99, the NS radius is assumed to be larger 
that 3i?s (where Rs = 2GM/c 2 is the Schwarzschild radius of a 
NS of mass M), but the accretion disc structure is not significantly 
changed up to the NS surface, and the radial velocity is always 
subsonic. If the radial velocity of accreting gas is supersonic at the 
surface (see e.g.[Pooham & NaravanlH99l for a possibility of the 
super sonic radial velocity in classical BL, and Kluzni ak & Wilsonl 
1991, for the "gap accretion" when the NS is within the innermost 
stable circular orbit), some fraction of the kinetic energy (associ- 
ated with a small radial velocity component) should be dissipated 
in an oblique shock, but most of it still remains stored in the kinetic 
energy of the gas rotating at the surfac e to be diss ipated later in the 
SL. The gap accretion model of Kluzniak & Wilson 1 1991) is rather 
similar to the SL, but they did not consider the fate of the accreting 
material and its spread over the surface. At very low accretion rate, 
both models should produce hard Comptonization spectra extend- 
ing to ~ 100 keV. 

The aim of this work is the calculation of the radiation spectra 
of the SLs and their comparison to the observed X-ray spectra of 
the BLs in the luminous LMXRBs. 



2 SPREADING LAYER MODEL 

The theory of the SL was developed by IS99 under the assumption 
of Newtonian gravity. They considered accretion of the pure hy- 
drogen plasma. Here we repeat the IS99 theory for plasmas of ar- 
bitrary chemical composition taking into account general relativity 
(GR) corrections using the pseudo-Newtonian potential. These cor- 
rections may be important, because the maximum effective temper- 
ature of the SL, which should be smaller than the local Eddington 
effective temperature TEdd, depends on the opacity and the gravity. 
The critical temperature is determined by the balance between the 
surface gravity and the radiative acceleration: 



GM 



(3) 



R?y/l-Rs/R c 

where a c — 0.2(1 + X) cm 2 g _1 is the electron scattering opacity, 
X is the hydrogen mass fraction, and <tsb is the Stefan-Boltzmann 
constant. It is clear, that solar (or larger) helium abundance together 
with the GR correction will lead to higher TEdd, and, therefore, to a 
higher maximum effective temperature of the SL. This could have 
important consequences for the determination of the neutron star 
parameters from observations. 

We use the pseudo-Newtonian potential in the form: 



*(r) = -c (l - Vl-ife/r) • 

This potential gives the correct GR surface gravity 

/ p\ GM 

flVi - Rs/R 

but gives the Keplerian velocity at the NS surface 



(4) 



(5) 




vi(R) = 



Figure 1. Geometry of the problem. 



GM 



(6) 



R^/l-Rs/R 
which is smaller than the correct GR value. 

2.1 Main equations 

Below we rewrite the IS99 equations for the SL for the pseudo- 
Newtonian potential J4j and considering arbitrary abundances. We 
consider the dynamics of the SL on the NS surface (see Fig. 0. 
The full hydrodynamic equation s, which describe this process are 
as follows (see for example lMihalasl 1 9781) . The continuity equation 
is 



+ V-(pv)=0, 



(7) 



dp 
dt 

where p is the plasma density, v is the vector of the gas velocity in 
the SL with components v v , vg and v r , which are velocities of the 
SL along longitude, latitude and radius correspondingly. Conserva- 
tion of momentum for each gas element is described by the vector 
Euler equation 



dv 



pv-Vv = -VP + f, 



(8) 



where P = P ra d + P s is the total pressure which is a sum of the 
radiation and gas pressures, and / is a force density. The energy 
equation for the gas in the SL is 



d_ 

dt 



1 2 

-pv + e 



1 2 

— pv 



f v-Vq + Q^ 



(9) 



Here e = e ra( j + e g is the total density of internal energy, where 
£rad = aT 4 is the radiation energy density and e g = (3/2)P g is 
the density of the internal gas energy. The first term on the right 
hand side is the power produces by the force density, the second is 
the energy, which is lost by radiation (q is a vector of the radiation 
flux), and the third is the heat, which is generated within a unit 
volume of the SL. 

Following IS99 we consider the steady state SL model in the 
spherical coordinate system (r, 9, if), where 9 is the latitude and if 
is the azimuthal angle (see Fig.0. We also assume that the SL has 
a small thickness (in comparison with the NS radius R, therefore 
the radial coordinate r — R), the radial velocity component is zero 



Spectra of neutron star spreading layers 3 



v r = 0, and it is axially symmetric (therefore, all of the deriva- 
tives 8/8<p equal to zero). In this case equations Q~j9j take the 
following form. The continuity equation is 

8 



1 



-(COS 8 pVg) = 0, 



Rcosddd 

the three components of the Euler equation are 

vl+vl \ __dP 



v e dvg v 



1 dP 
R~89 



+ fe 



ve 



8v^ 



and the energy equation is 



1 



8 



R cos 8 86 



cos 8 ve + e + P 

8q 

feve + f v v v - — + i 
or 



(10) 

(11) 
(12) 
(13) 

(14) 



where 



2 2.2 

«o = v v + v e ■ 



(15) 



Here the radiation flux has only one (radial) non-zero component 
and its divergence is computed in the plane-parallel approximation 
which is a consequence of our assumption of small height of the 
SL. A small azimuthal component of the radiation flux arises due 
to the aberration, which we neglect here. 

It is clear that equation <1 H can be solved independently on 
equations I12l - il3| and we can consider some averaging over the 
layer's height. Therefore, we arrive at a one-dimensional problem. 
In this case instead of equations 1101 — 1141 we obtain 



8 



R cos 6 89 



cos 8 



J pvedrj = 0, 



(16) 



Pdr+R / /odr,(17) 



f 8v<p 

J pV0 ^e dr 



tan# / pv<pVedr = R / f^dr 







R cos 8 86 



ve | -pv + e 



f g v g dr+ / f v v v dr 



P dr 



(18) 



(19) 



h dr. 



Here the integration over radius is from R to R + hs , where hs is 
the local SL thickness. We define the corresponding force densities 
in the next section. 



2.2 Vertical averaging 

The one-dimensional equations for the SL structure are derived us- 
ing the averaging along the height at a given latitude. It means 
that we have to calculate all the integrals in equations <16t - <19| 
for some model of the SL vertical structure. The simplest way is 
just to consider the variables averaged over the height. 

IS99 used a more complicated model of averaging. They con- 
structed a simple model of the SL using assumptions that veloc- 
ities ve, v^, and the radiation flux do not depend on the height 
q{r) = const = crsBTeff • The latter suggestion means that all of 
the thermal energy in the SL is generated at the bottom. This model 



is described by the hydrostatic equilibrium equation i 1 1 i taken in 
the form 



dP 



dm ^ = 9° R ' 
and the radiation transfer equation in the diffusion approximation 

de rad 3q 



2 i 2 
V^+Vg 



(20) 



dm 



(21) 



Here and below we use a new independent variable: a column mass 
m and a new geometrical coordinate z, which are defined as 



dm = pdz = —pdr. 



(22) 



Coordinate z has an opposite direction relative to r and z=0 at r — 
R + hs . We also defined the r component of the force density 

fr = -gop- (23) 

Equations I20t — 12 1 i have to be supplemented by the equation of 
state 



_ pkT e ra d 
/im p 6 



(24) 



where p, — 4/(3 + 5X) is the mean molecular weight and m p is 
the proton mass. 

Equations 1201 — 1241 can be solved with the simple boundary 
conditions P(m = 0) = 0, T(m = 0) = 0: 



g e sm, 

3q 



Pcor i CtC 3 

p — pm p — — I m 



k V 3qa c 



1/4 



T = I —ma. 



1/4 



1/4 



(25) 
(26) 

(27) 
(28) 



where r c = ma c is the electron scattering optical depth of the 
layer, and 



<7cor = ffcff — ffrad = ffefi 0" e . 

C 



The column density m is related to the geometrical depth z 

_ (pm p g, ' 



cor I (XC 4 



(29) 



(30) 



4 4 a c fc 4 3q 

which gives the following dependence of temperature on height 

(31) 



= pm p g cov 



4A: 

Following IS99, we consider the values of temperature and 
density at the bottom of the SL Ts and ps as parameters. In this 
case the local SL thickness hs is: 

4fcT s 



hs 



(32) 



pm p g COI 

We can also calculate all of the integrals in equations <16M19> : the 
total surface density 



pdz = m(z = hs) = Eg, 
the pressure surface density 

q Pdz = iffeff/isSs = \ gcff SsfcTs, 
5 5 pm p g COT 



(33) 



(34) 
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the surface density of the internal energy 



M = 4nR cos 9 v g Y, s 



(40) 



the local flux 



£rad + ~Pg I dz 



cff + flrad 



Es^s 



^4 ac r s 4 
q = asBJcff = ^ — =r, 
ocr c 2jS 

and the enthalpy flux 

H — -Eint + 



Pdz 



5 , 3 \ E s fe 



(35) 



(36) 



(37) 



If we take X — 1, all these relations will be the same, as derived 
by IS99 with one exception: there is no potential energy of the SL 
in our energy equation. Thus our expression for H contains a factor 
5/2 instead of 7/2 as in IS99. Below we will show that this pro- 
duces only a small quantitative differences between our and IS99 
models. 

In the IS99 model there are two forces, which give contribu- 
tion to the force density in equations <16M19t . These are the grav- 
ity force, which has only the radial component (see above) and the 
force arising due to the friction between the NS surface and the 
SL. This force is directed along the NS surface and is expressed in 
the IS99 model through stress r and its azimuthal and meridional 
components t v and rg. IS99 have parameterized it in the form: 



TO 



f v dz = a h psv<pVo, 



fedz = a b psvevo, 



(38) 



where Qb = w* /vq is the parameter of the stress parametrization, 
w» is the velocity of turbulent pulsations. We should note that a?b 
is not the same a that is used in the accretion disc theory. In accre- 
tion discs, a (in the first approximation) is the square of the ratio 
of the turbulent velocity to the sound speed a = wj/ef and can be 
quite high, up to 0. 1-1 . In the SL, the plasma velocity vq is close to 
the Keplerian velocity at the NS surface and is orders of magnitude 
larger than the sound speed. The velocity of turbulent pulsation is 
also limited by the radiation viscosity at the SL bottom. IS99 care- 



fully investigated this matter and estimated Qb 



10" 



We used 



this value in most of the paper. 

IS99 have ignored the mechanical work between the SL and 
the NS (which accelerates the stellar rotation). In our work we use 
the same approximation. A fraction of the kinetic energy of the ac- 
creting gas that goes to increase the rotational energy of the NS is 
approximately 2£1,/£Ik, where fi„ is the NS angular velocity and 
Ok is the Keplerian angular velocity at the NS surface. As we con- 
sider a non-rotating NS and the characteristic time to increase its 
angular velocity is orders of magnitude larger than the character- 
istic time of the SL t = R/v e = 10 km/10 3 km s" 1 = 0.01 s, 
ignoring the mechanical work on the NS is a reasonable approxi- 
mation. In this approximation, therefore, all the work due to friction 
between the SL and the NS transforms to heat: 



dz — ~(t v v^, + rgvg) — tvq. 



(39) 



2.3 One-dimensional model of the spreading layer 

Using relations <32> - <39t we can rewrite equations, which describe 
the one-dimensional SL structure. The continuity equation can be 
rewritten via the accretion rate as 



Therefore, the product cos 6 vgT,s = const. The Euler equations 
arc 



a v ' i 4 a ( 9e« kTs „ 

COS 9 Vgl^sVg + — COS 6 E 



5 \ #cor M m P 

~Rcos9 rg — sin# Esv^,, 
Es«e(cos Dtp) — —Rcos8t v . 



(41) 



(42) 



Here the prime means the derivative over 9. The second term in the 
left hand side of equation J4 1 1 is the lateral force gradient, and the 
terms in the right hand side are the components of the stress force 
and the centrifugal force. We see from equation 1421 that the mo- 
mentum along if coordinate is changed due to the friction with the 
NS surface only. The system of equations is closed by the energy 
equation 



Es^e 



Wo , 2 kT s 5g cB + 3g rad 



2^5 urrip 



ffcor 



-Rq. 



(43) 



The system of equations 1401 — 143^ can be transformed to the 
three dimensionless equations for v lp (8), vg{9) and Ts(9) as was 
done by IS99. These equations are solved with the boundary condi- 
tions at the transition zone between the accretion disc and the SL: 
the initial latitude, where the SL starts, is close to the NS equator 
#o ~ 10 -2 ; the initial relative deviation S of v v from the Kep- 
lerian velocity v v — hk(1 — 5); and the initial ratio of the ki- 
netic energy of the SL along 9 coordinate and it's thermal energy 
9 ee (um p vj )/kT s . 

As was demonstrated by IS99, the solution of equations 1401 — 
(35) depends very little on 9o (if #o sufficiently small < 0.1, see 
below) and 5, but strongly depend on parameter Q. We choose 
the solutions which are closest to the critical solution (in this so- 
lution vg is equal to the sound speed at the maximum latitude of 
the spreading layer), but slightly subsonic. The necessary critical 
value of parameter O is found by the bisection method. 

The distributions of v$, the effective temperature T e g, and the 
surface density Es along the latitude 9 for four models with the 
same accretion rate, corresponding to first model luminosity are 
shown in Fig. [2] The first model (shown by the solid curves) is 
our model with the pseudo-Newtonian potential and solar hydro- 
gen abundance X = 0.7. The dashed curves are for our model 
with GR correction, but with pure hydrogen X = 1; the dotted 
curves correspond to our model without GR corrections (Rs = 
in the equations) with solar hydrogen abundance (X = 0.7), while 
the dot-dashed curves are for the IS99 model with GR corrections 
and solar hydrogen abundance. It is clear that the solar abundance 
lead to a narrower spreading layer with a smaller surface density. 
A higher helium abundance as well as the the GR corrections lead 
to a higher effective temperature and a larger latitudinal velocity. 
Our model gives a slightly wider SL with slightly smaller latitu- 
dinal velocity but same effective temperature and surface density. 
Most calculations below were performed for our model with the GR 
correction and solar abundances. The surface density and the effec- 
tive temperature distributions along the latitude for models with 
0.1, 0.2, 0.4 and 0.8 of the Eddington luminosity are presented in 
Fig. [5] Variations of parameter Qb lead to some changes in the SL 
structure. The SL column density is inversely proportional to a b , 
while the resulting effective temperature decreases only by about 1 
per cent with decreasing of Qb by an order of magnitude. 

The lower boundary of the SL was taken very close to the 
equator, 9o « 0.01 — 0.001, in IS99. Formally, the accretion disc 
thickness is close to zero at the inner boundary, if we take the usual 
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Figure 2. The distributions of (a) vg, (b) the surface density £3 and 
(c) the effective temperature T e g along latitude 9 for four models. Solid 
and dashed curves are for our model corrected for the GR using pseudo- 
Newtonian potential and X = 0.7 and X = 1, respectively. The dotted 
curves is our model for the Newtonian gravity and X = 0.7, while the dot- 
dashed curves are for the IS99 model corrected to GR and X = 0.7. In all 
cases M = 1.4Mq and R = 12 km. The accretion rate of all models is the 
same, corresponding to 0.2L Edd of the first model. 



inner boundary condition for the component of the stress tensor 
W r <p{Rin) = 0. In the case of the accretion disc around a NS this 
condition is not correct, and the disc thickness at the NS surface 
is considerable. The luminous accretion disc half-thickness can be 
evaluated from the balance of the radiation force and ^-component 
of gravity: 

3(7, 



8ttc 



-M. 



We calculated the SL models with the initial latitudes 9oi = 
arcsin(zo/-R) and #02 = arcsin(0.5;Zo/-R). The surface density 
and the effective temperature distributions along the latitude for the 
second case are shown in Fig.|4]. The qualitative behavior of these 
distributions is close to the case of small 9o with some shift along 
the latitude. The main difference is the maximum possible luminos- 
ity of the SL. At this luminosity the SL reaches the NS poles. For 
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Figure 3. The distributions of the surface density Ss and the effective 
temperature T c g over the latitude 9 for the SL models with different lumi- 
nosities. Initial SL latitude is 9q = 10 — 2 . 

#02, the maximum possible luminosity is about 0.4 L^dd, while for 
f?oi it is about 0.25 LEdd- 



2.4 Vertical structure of the spreading layer 

The IS99 SL model was constructed under an assumption that the 
local layer velocity and the radiation flux along height is constant. 
It means that the SL is decelerated and energy is liberated in the 
infinitely thin layer at the NS surface. This is an approximation 
only and the velocity distribution should not be uniform and the 
energy should be generated at all heights. Therefore, we need a 
more detailed vertical model of the SL for calculating its radiation 
spectrum. 

For evaluation of the viscosity parameter a^, IS99 used classi- 
cal theory of hydrodynamic BLs lLandau & LifshitJl959h with the 
logarithmic velocity and the energy generation distribution along 
the height. In this case, both the energy generation rate and the ve- 
locity gradient are inversely proportional to the distance from the 
NS surface z 

JL K oc _. (45) 

dz dz z 



(44) oc ^ 



It is clear that these dependencies cannot be correct in a SL. The 
SL has a finite thickness with a low density at the surface. How- 
ever according to equation 1451 some amount of energy has to be 
generated in the surface layers. 

At present time, a theory of the radiation-dominated turbulent 
boundary layer does not exist. Thus we here can only make simi- 
lar assumptions about the energy generation and velocity gradient 
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Figure 4. The distributions of the surface density Es and the effective 
temperature T c g along latitude 9 for SL models with different luminosities. 
Initial SL latitude 8q = arcsin(0.5zo/-R). 



along the height. We assume that these values are inversely propor- 
tional to the surface density measured from the NS surface: 



dm Es 
dv 

— A 

dm Es 



go 



'(.'(> 



(46) 
(47) 



where A = 2.5a^f 2 , qo and vrj are the local radiation flux and the 
average layer velocity at a given latitude obtained from the one- 
dimensional model. Equations I46i and 1471 are very close to the 
IS99 SL model assumptions (the layer is decelerated and the en- 
ergy is generated at the bottom of the layer). Integration of these 
equations yields 



q(m) = q Q 



I + 1 In ( 1 — — 



(48) 



v(m) 



Vo 



1 + A In 



Es 



(49) 



These equations can be used up to some critical column density 



m, = E s (l-exp[-4 ']) , 

which is very close to the local surface density. 

The hydrostatic equilibrium equation then reads 



dm 



90 



v 2 {m) q(m) 



R c 
and the radiation transfer equation is 



(50) 



(51) 



1 de q(m) 



(52) 



3 dm c 

The temperature and the gas pressure distributions along the height 
are: 



T(m) = T ci 



mCTe [ 1 — A 



1 + 



E s - m 



(53) 



'"""El 



1/4 



-£m 1 -2A + 2A 1 
R 



(54) 



A In [ 1 - — ) - 2A + 2 



c 



in - A [ m + (Es - m) In ( 1 - 



These solutions are obtained using the boundary conditions at the 
surface e(m = 0) = 2qo/c and P s (m = 0) = 0. 

At the same time, there is a disagreement between this verti- 
cally explicit model and one-dimensional model, because the ve- 
locity and the flux vertical profiles are different. We suggest, that 
the model can be made self-consistent, if we find a new value of the 
surface density Eg at a given latitude, which conserves the mass 
flux 



«oEs = 



si 



v(m)dm = «q(1 — j4)E s . 



(55) 



Therefore, the new value of the surface density E s = Eg/(1 — A). 
For «b = 10~ 3 , this gives Eg = 1.086Es and we take these val- 
ues below for our calculations. There are similar disagreements for 
other integrals over the height in equations 1 16i — 1 1 91 . For example: 

^oEs = / S w 2 (m)dm = wo(l-2^ + ,4 2 )Es. (56) 
Jo 

In this case, we have to take a new value of the surface den- 
sity E'g = E s /(1 - 2A + A 2 ), which gives E' s = 1.171E S if 
Qb = 10~ s . Our vertically explicit models disagree with the one- 
dimensional ones by about 10 per cent. Fortunately, the emitted 
local spectra depend very little on the surface density of the SL. 



3 SPECTRUM OF THE SPREADING LAYER 

For calculation of the SL spectra we divide it into a number of rings 
over the latitude which have different effective temperatures T e g, 
matter velocities Vo and surface densities Es. We then calculate the 
vertically explicit model for each ring, solve the radiative transfer 
equation and obtain the local SL spectrum. Then we integrate local 
spectra from the SL surface accounting for the general and special 
relativity effects. 



3.1 Local spectra 

To calculate a vertically explicit hydrodynamical model with the 
radiation transfer w e use s tandard methods for stellar atmospheres 
modelling iMihalas 1978). Our models are obtained in the hydro- 
static and the plane-parallel approximations. The effective temper- 
atures of the considered SL models are rather high (~ 2 keV) and 
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these models are similar to the atmospheres of bursting NSs, where 
Compton scattering have to be taken into account. 

The vertically explicit local SL model is described by the fol- 
lowing equations: the equation of hydrostatic equilibrium 15 1 i . the 
energy generation law I48i . the velo city law 1491. the RTE account- 
ing for the Compton effect using the Komcaneets] ll957l) operator: 



dx 



dJ v 
dx 



(J„ - B v ) - 



kT 



k v + £T e 



or, T E fi 



1 + 



m c c 
CJ V 



(57) 



where x = hv/kT c g is dimensionless frequency, f v (r v ) ~ 1/3 is 
the variable Eddington factor, J v is the mean intensity of radiation, 
B v is the black body (Planck) intensity, k v is the opacity due to the 
free-free and bound-free transitions, a B is the electron (Thomson) 
opacity, T is the local electron temperature, T e g is the effective 
temperature of SL at a given latitude, and C = c 2 h 2 /2(kT c g) 3 . 
The optical depth r v is defined as 



(k y + <r c )dm. 



(58) 



These equations have to be completed by the energy balance equa- 
tion 

kT 



k v (J v - B u ) dv - -^—z^- - cr e 

4tt dm m c c z 



X 



(59) 



J v du — 



Teg 

T 



< /, ! 1 + < —^- ) dv 



= 



and by the ideal gas law 

P<? = iV tot fcT, 



(60) 



where JV to t is the number density of all particles, as well as by 
the particle and charge conservation laws. We assume local ther- 
modynamical equilibrium (LTE) in our calculations, so the number 
densities of all ionization and excitation states of all elements have 
been calculated using Boltzmann and Saha equations. 

For solving these equations and computing th e local SL model 
we used the Kurucz's code ATLAS lKuruczll970l,l 199*31) modified 
for high temperature. All ionization states of the 15 most abundant 
elements are taken into consideration. The photoionization cross- 
sections from the ground states of all ions are calculated using PH- 
HT 2 code jVerner et alll99d) . For details see lSwartz et all 1 2002) 
and Ib ragimov et alJl20 03 ). The code was also modified to account 
for Compton scattering. 

The scheme of calculation is the following. First, the input 
parameters of the local SL model are defined from the total one- 
dimensional SL model (see Sect. 12.31 : the effective temperature 
TeB, the surface gravity go, the surface density Es, and the local 
average layer velocity Wo. Then the analytical vertically explicit 
model I53H55I are calculated together with the new value of sur- 
face density E^ = Eg/(1 — A). The calculations are performed 
for the set of 98 column densities m, distributed logarithmically 
with equal steps from m = 10~ 5 g cm~ 2 to 0.99m, . The gas pres- 
sure, which is found from equation I54i . is not varied during the 
iterations. 

For this starting model, all number densities and the opacities 
at all depth points and all the frequencies (we use 300 logarithmi- 
cally equidistant frequency point s) are calculated. The RTE 1571 is 
solved by the Feautrier method (Mihalas 1978; Zavlin & Shibanov 
Il99lt iPavlov et alJll99ll: iGrebenev & Sunvaev|l2002h iterativelv. 
because it is non-linear. Between the iterations we calculate the 
variable Eddington factors /„ and h v , using the formal solution of 



the RTE for three angles. Usually 5-6 iterations are sufficient to 
achieve convergence. 

We used the usual condition at the outer boundary 



d(f v J v ) 



(61) 



(62) 



where hv is the surface variable Eddington factor, and the inner 
boundary condition 

dJ v _ dB v 

The outer boundary condition is found from the lack of the in- 
coming radiation at the SL surface, and the inner boundary con- 
dition is obtained from the diffusion approximation J v ps B v and 
q v ~ 4tt/3 x dB v /dr v . This condition is satisfied for any SL 
optical thickness, because the SL bottom is the NS surface. The 
boundary conditions along the frequency axis are 

Jv = B v (63) 



at the lower frequency boundary, v = u m i n 
0.03 eV < feTefr) and 

dJv „ r T cS I CJv\ 
X ^x--^ + — Xj »\ 1 + — ]=Q 



10 14 Hz (hl> m in 



(64) 



at the higher frequency boundary v = f max = 3 10 19 Hz 
(ta ma x w 100 keV » feTeff). Condition <63) means that at the 
lowest energies the true opacity dominates the scattering k v 3> c e , 
and therefore J v w B v . Condition <64> means that there is no pho- 
ton flux along the frequency axis at the highest energy. 

The solution of the RTE <57t should also satisfy the energy 
balance equation 1591 and the surface flux condition 



q v {m — Q)dv — assT eB . 



We calculated the relative flux error along the depth 



(65) 



q(m) 



e q (m) = 1- r^","'., , (66) 
Jo qv(m)dv 

where q(m) is found from the energy generation law <48> . and 
q v {m) is radiation flux at a given depth obtained from the first mo- 
ment of the RTE 

4tt — = qv. (67) 

Then the temperature corrections were evaluated using three dif- 
ferent procedures. The first procedure is the integral A-iteration 
method based on the energy balance equation 1591 which was mod- 
ified to account for Compton scattering. It works well in the upper 
layers. The second one is the modified Avrett-Krook flux correc- 
tion, which uses the relative flux error and is good in deep layers. 
And the third one is th e surface correc tion, which is based on the 
emergent flux error. See lKumczlil97bT) for the detailed description 
of the methods. 

The iteration procedure is repeated until the relative flux error 
is smaller than 1 per cent, and the relative flux derivative error is 
smaller than 0.01 per cent. As a result we obtain the self-consistent 
local SL model together with the emergent spectrum of radiation. 

Our method of calculation was checked on the atmosphere 
model of bursting NS. The equations which describe the bursting 
atmosphere are simpler, because there is no velocity field along 
the surface (no = 0) and the integral flux is constant along depth 
(dq/dm — 0). We compared our model atmospheres with the 
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Figure 5. (a) The spectra of the bursting NS model atmospheres with the 
effective temperature T e g = 2 10 7 K, solar hydrogen/helium abundances 
without heavy elements and log g = 14.2 and 15. The black body spectrum 
corresponding to the effective temperature is shown by the dotted curve, 
(b) The corresponding temperature structure as a function of the Thomson 
optical depth. The dashed curves show the spectrum and the temperature 
structure for one model computed not accounting for Compton scattering. 



most recent models of iMadei. Joss & Rozanskal l2004t) . The radi- 
ation spectra and the temperature structure for some models with 
Teg = 2 10 7 K, solar H/He abundances, and various surface grav- 
ities are shown in Fig. [5] T hese result s are in a perfect agreement 
with the results of Madei et al. (2004). The emergent spectra and 
the temperature structure for the models with the solar abundance 
of heavy elements are shown in Fig. [6] 

In the surface layers, local cooling is small because of the low 
density, and the temperature equals the Compton temperature of 
radiation which is slightly higher than the effective temperature. In 
dipper layers, at r c ~ 0.1, the cooling due to thermal emission 
(free-free and bound-free) becomes important (as thermal emissiv- 
ity per gram is proportional to density) and the temperature de- 
creases. At large optical depth the temperature rises again and fol- 
lows the Te relation, typical for a grey atmosphere. At higher sur- 
face gravity (at fixed T e g ), the plasma density is higher, resulting 
in a more significant temperature dip. We see that heavy elements 
have rather minimal influence on the models close to the Eddington 
limit (lower g). 

The comparison between the bursting NS models and different 
local SL models for the same T e s and effective log g is shown in 



Figure 6. Dependence of the bursting NS model atmospheres on metal 
abundance. Dashed curves correspond to the solar abundance of metals and 
solid curves to the zero metal abundance. Examples are for solar H/He abun- 
dances, the same effective temperature T e g = 2 10 7 K and two values of 
log g = 14.2 and 15. Black body spectrum corresponding to the effective 
temperature is shown by dotted curve. 



Fig.0 F°r this SL model we use the vertical structure model, which 
is described in Section l2~4l We also investigated, whether the model 
for the vertical structure is important for the emergent spectra of 
local SL models. We also calculated the SL model with the constant 
velocity and flux derivatives: 

dv vo 
dm Es 
and 



(68) 



(69) 



dq _ go 
dm Es 

In this case the mass flux conservation requirement <55t leads to 
Eg = 2Es. The spectrum and the temperature structure of this 
model are shown in Fig. Q by squares. The surface temperatures 
of the local SL models are higher than the bursting NS model sur- 
face temperature. The reason is the non-zero flux derivative in the 
energy conservation equation I59i . This means that a part of the 
energy is released in the upper atmosphere and is heating it ad- 
ditionally. The smaller the surface density (i.e. the larger the flux 
derivative), the higher the surface temperature. But the differences 
in the temperature structure have very small influence on the emer- 
gent spectra. Therefore we conclude, that details of the vertical 
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Figure 7. Comparison of the bursting NS model (dashed curves) and dif- 
ferent local SL models for the same T c fj and effective gravity log g. Solid 
curves correspond to the local SL model (with the vertical structure model 
described in section 12.41 with surface density Eg = 630 gem -2 . The 
spectrum and the temperature structure of the local SL model with the con- 
stant velocity and flux derivatives are shown by squares. Dotted curves cor- 
respond to local SL model (Sect. 12.41 with Eg = 63 g cm -2 . 

structure have negligible influence on the emergent spectrum for 
the optically thick models (Ss > 100 g cm -2 ). 

It is well known that the model spectra of bursting NS close 
to the Eddington limit are well described by a diluted Planck spec- 
trum with the color temperature T c = f c T e g with the hardness 
factor f c varyin g in the inter val 1.6-1.9 and the dilution factor 
D = /T 4 . |PavloveTal]h99ll) have derived an analytical formula 
for the hardness factor, which successfully describes high luminos- 
ity (L « -Lridd) burst spectra: 

f c = (0.15 InCi + 0.59)~ 4/5 C* /ls e 3/20 , (70) 

where Ci = (3 + 5X)/(1 - t) and £ = L/L E dd = ffrad/ffeff- 
Equation 1701 works well also for models with heavy elements. 

The local spectra of the optically thick SL (with L > 
0.2Z/Edd) are very similar to the burst spectra with corresponding 
parameters (see Fig. |7J. The local SL are very close to the Edding- 
ton limit 

5rad = — CT e « ffO - -5- = fleff • (71) 

c R 

For example, the distributions of the ratios g r ad/<?o and g Ta ,d/g B s 
along the latitude for SL models with three different luminosities 
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Figure 8. (a) Radiative and effective accelerations and (b) the hardness fac- 
tor f c (b) as functions of latitude for the SL models of different luminosities 
L/ L Edd . 



are shown in Fig.[8^. Corresponding hardness factor distributions 
are shown in Fig. [8j>. The comparison of the two local SL spec- 
tra (close to equator and at higher latitude) with the diluted Planck 
spectra and hardness factors given by equation I70i are shown in 
Fig.|5^. Closer to the equator, effective gravity is low as centrifugal 
force is large. The gas is levitating above the NS and £ is close to 
unity. The energy dissipation and the effective temperature are low. 
Thus, fc is large and the spectrum is close to the diluted Planck. 
At higher latitude, the layer is decelerated, while the energy dissi- 
pation and Tog grow. However, the effective gravity grows faster 
reducing £ and the color correction f c . The spectrum shows devi- 
ations from the diluted Planck spectrum at low energies. At high 
energies, the Wien part of both spectra can well be described by the 
diluted Planck. 



3.2 Integral spectra 

Now we can compute the integral total model spectrum of the SL, 
which is seen by a distant observer accounting for the relativis- 
tic effects such as gravitational redshift, light bending, relativis- 
tic Doppler effect and aberration. We take into account only half 
of the SL because another half is hidden by the accretion disc 
and divide the SL surface on 10 latitude rings and on 100 angles 
in azimuth. In a spherical coordinate system, where the accretion 
disc coincides with the 8 = 0° plane, the spectrum of the SL is 
IPoutanen & Gierlinskl2003t) 
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Figure 9. (a) Local SL model spectra with high and low effective tempera- 
tures (solid curves) compared to the diluted Planck spectra with f c given by 
equation 1701 . (b) Integral spectra (isotropic luminosities) of the same SL 
model (solid curves) for different inclination angles i = 0° (softer spec- 
trum) and 90° (harder spectrum) compared to the total spectra obtained by 
integration of the diluted Planck spectra (dashed curves). 

Fe = Jyi J J V 3 S a I(E', cos a, 9) cos a cos 9d9 dip. (72) 



Here the observed and the emitted photon energies are con- 
nected by the relation E = E' T) 5, where 77 = \/l — Rs /R, 
the Doppler factor 8 = 1/7(1 — /3cos£), f3 = v v (9)/c (here 
we neglected low latitudinal velocity), the Lorentz factor 7 = 
l/\/l — P' 2 , and cos£ = —sin a sin i sin ip/ sin ib. Th e light 
bending is accounted for by the relation jBeloborodovl2002l) 

cos a = + rj 2 cos ip, (73) 
R 

where cos ij) = cos i sin 6 + sin i c os 9 cos ip, and the rela tivistic 
aberration gives cos a' = 5 cos a JPoutanen & Gierlinskill20()3T) . 
Here i is the inclination angle of the NS polar axis to the line of 
sight, D is the distance to the observer, and #sl is the SL boundary. 
Only visible surface elements with cos a > give contribution to 
the total spectrum. 

The emitted specific intensity I (E' , cos a' , 9) is taken from 
the computed local SL flux assuming angular dependence for the 
electron scattering atmosphere 

I(E',cosa',9) = qw ^ (0.4215 + 0.86775 cos a) . (74) 

This formula gives a good approximation to the specific intensity 
of the emergent radiation (see Fig. llOi . 
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Figure 10. Specific intensities of the emerging radiation of the local SL 
model (here [i = cos a') with the same parameters as in Fig.Qcomputed 
exactly (solid curves) and using formula [TU. 
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Figure 11. Dependence of the SL integral spectra on the inclination angle 
i to the line of sight (a) and on the luminosity (b). 



The total spectra of the SL model for two inclination angles, 
i = 0° and 90°, are shown in Fig.|9j). The spectra computed using 
the local diluted Planck spectra are shown also for comparison. The 
difference is very small in the high energy part (E > 10 keV) and 
more significant at lower energies (E < 7 keV). 

Dependence of spectral shape on the inclination angle is not 
significant (see also Fig. II lh ). Differences between the SL spectra, 
which are seen at different inclinations are comparable to the dif- 
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Figure 12. Comparison of the observed spectra of the BLs of LMXRBs 
obtained by the Fourier-frequency resolved spectroscopy and the model 
spectra of the SL at two inclination angles to the line of sight. Filled cir- 
cles give the spectrum of GX 340+0 in the normal branch (see Fig. 12 in 
iGilfanov et al . 2003 ), op en circles correspond to the spectra of five Z- and 
atoll sources (from Revnivtsev & Gilfanov 2006). 



ferences due to change in the SL luminosities (see Fig. ll lb ). It is 
interesting, that the total spectra can also be well described by the 
diluted Planck spectrum. 

The color temperature depends slightly on the assumed tur- 
bulence parameter ctb. Decreasing «b by an order of magnitude 
increases T c by 0. 1 keV. 



4 COMPARISON WITH OBSERVATIONS 

In LMXRBs a weakly magnetized NS is surrounded by the accre- 
tion disc which transforms to the boundary/SL close to the NS sur- 
face. At present, about 100 LMXRBs are known. They can be di- 
vided into two different classes. The Z-sources are very luminous 
(L ~ 0.1 — 11/Edd) and have relatively soft, two-component spec- 
tra. Both components are close to the black body with color tem- 
peratures of about 1 keV and 2-2.5 keV. The atoll sources are less 
luminous (L ~ 0.01 — 0.05I/Edd) and are observed in two states, 
the high/soft and the low/hard. In the soft state, the radiation spectra 
are similar to those of the Z-sources, while in the hard state they are 
close to the spectra of the Galactic black hole sources in the har d 
states (e.g. Cyg X-l, see e.g. IPoutanenll 19981: iBarret et al-fe oOO). 
These hard spectra are well described by unsaturated Comptoniza- 
tion of soft photons in the hot (kT ~ 30 - 100 keV) optically thin 
(r e ~ 1) plasma. 

The soft component can be associated with the radiation 
from th e accretion disc, while the hard one with the boundary/SL 
iMitsuda et all 19841) or possibly with a corona or hot optically thin 
inner accretion flow (see discussion in Done & Gierlihski 2003) 
in case of low-luminosity atoll sources. At high luminosities, the 
BL is optically thick and its effective temperature is higher than 
that of the accretion disc, because the BL is smaller than the ac- 
cretion disc, while their luminosities are comparable. The hard 
component is also more variable than the soft componen t at the 
timescales from millisecond to 1000 s econds i Mitsuda et all 1 98*41 
Gilfa nov". Revnivtsev & Molkovlbooll) . The Fourier-frequency re- 
solved spectroscopy confirms that a component variable at high 
frequencies (and sometimes showing quasi-periodic oscillations, 



sec lvanderKlisll200d) has a blackbody-like spect rum with the 
color temperature T c = 2 .4 ± 0.1 keV iGilfanov et all l2003t 
iRevnivtsev & Gi lfanov 2006) which is very similar for the five in- 
vestigated sources. On the other hand, the variability of the soft 
component is very similar to the variability of the soft component 
of black hole sources in their soft states, which is associated with 
the accretion disc. Based on these arguments, we associate the hard 
blackbody-like component with the BL and compare our theoreti- 
cal SL spectra with it. 

Spectra computed for one SL model together with the ob- 
served BL spectra obtained by the Fourier-frequency resolved spec- 
troscopy ( Gilfan ov et al J 120031 IRevnivtsev & Gi lfanov 200^) are 
shown in Fig. 1121 We see a very good agreement between theoreti- 
cal spectra and the spectrum of GX 340+0 at the normal brach (at 
high accretion rates). The spectra of five Z- and atoll-sources (open 
circles) are similar to our SL spectra at high energies, but have a 
soft excess. This excess may be related to the emission of the clas- 
sical BL, the inner part of the accretion disc. The observed spectral 
similarity gives us a confidence to try to determine NS parameters 
from the observed spectra. 

As we have shown above the spectrum of the SL can be rep- 
resented by a diluted blackbody. The effective temperature of ra- 
diation is determined by the critical temperature from equation 

where the left-hand side is multiplied by £, the ratio of the 
local flux to the critical Eddington one (reduced due to the ac- 
tion of the centrifugal force). The observed color temperature is 
T c = f c \J 1 — Rg I R TEdd, where corrections are made for spec- 
tral hardening and gravitational redshift. For the known color cor- 
rection and £, the NS radius as a function of compactness M/R can 
then be found from 
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1- — 
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(75) 



2a SB TcV e R 

Assum ing f c — 1.6 — 1.8 and I — 0.8. IRevnivtsev & Gilfanovl 
(2006) obtained constraints on the NS mass-radius relation (shown 
in Fig. ll3l bv dotted curves). The maximum NS radius is reached 
forR s /R = 2/5: 
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(76) 



Here instead we calculate exactly a grid of the SL model spec- 
tra, where the main input parameters are the NS mass M and radius 
R, and the SL luminosity. The NS mass is varied from 1 to 2 Mo 
with a step of 0.2 Mq, and the NS radius is varied from 10 to 24 
km with a step 1 km. Only the models with R > 3Rs are consid- 
ered. We take Qfb = 10~ 3 and luminosity of 0.4LEdd, and compute 
spectra for four inclination angles 0, 30, 60 and 90 degrees and for 
three chemical compositions: pure hydrogen (X — 1), solar abun- 
dance (X = 0.7), and pure helium (X = 0, Y = 1). The spectra 
are fitted by the black body and the corresponding color tempera- 
tures are found. 

Models with higher He abundance have a smaller hardness 
factor as can be seen from equation I70i . However, the local effec- 
tive temperature of the SL is higher for larger He abundance (see 
eq.|5|and Fig.[2j;). The higher T c fr leads to a higher color temper- 
ature of the integral SL spectrum. For example, at NS radius of 13 
km and mass 1.4M0 pure hydrogen models give color temperature 
of about 2.5 keV, while pure helium models produce harder spectra 
with T c w 3 keV. 

Contours corresponding to the color temperature equal 2.3 
(right), 2.4 (central) and 2.5 keV (left) are shown on the Al — R 
plane (Fig. ll3> together with the NS models for various equations 
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Figure 13. Permitted region (shaded) for the NS 
mass and radii, which can have SLs with color 
temperatures 2.4 ± 1 keV iGilfanov et"aT] i2003; 
Revnivtsev & Gilfanov 2006) at luminosities similar 
to those observed in the Z-sources. The left shaded 
region (with boundaries and the centre marked by 
circles) is for the accreting matter containing hy- 
drogen only, the middle one (boundaries and centre 
marked by squares) is for the solar composition, and 
the right one (marked by triangles) is for pure he- 
lium. The left boundary of each region corresponds 
to T c = 2.5 keV, while the right boundary gives 
2.3 keV. Dotted curves correspond to a simple estima- 
tion of M — R relation given by equation 1751 with 
f c = 1.7, £ = 0.8 for the three chemical compo- 
sitions. Various theoretical mass-radius relations for 
neutron and strange stars are shown f or comparison 
lHaensel. Potekhin & Yakovlevl UooS) . Dotted hori- 
zontal line corresponds to the NS mass of 1.4Mq. 
The almost vertical dashed curve corresponds to the 
thermally emitting NS with the apparent radius of 
Roo = 16.5 km (see eq. 1771 and I riiiiii>cl l2005ll . 
The lower limit on the NS mass as a function of radius 
for a given rotational period P (in ms) as derived from 
the NS stability M/My > 0.865(i?/10 km) 3 P^ s 
lLattimer & Prakashl200l) is presented for P = 1.6 



of state. These iso-temperature curves are shown for the inclination 
angle i — 45°. Comparison of the observed spectra to the theoreti- 
cal spectra of the SL constrains the NS radius at 13.5 ±1.5 km (for 
pure hydrogen X — 1 model), 14.8 ± 1.5 km (solar composition 
X = 0.7) and 19 ± 1.5 km (pure helium X = 0, Y = 1) assuming 
the NS mass of 1.4 solar mass. For pure hydrogen and solar abun- 
dance, the permitted radii are consistent with the hard equation of 
state of the NS matter. If the composition is solar, but the heavier 
elements are able to sink, the emitted spectra would correspond to 
a pure hydrogen atmosphere requiring thus smaller radii. 

Increasing the inclination to 90° increases the deduced NS 
radii by about 10 per cent, while assuming i = 0°, gives a 15 per 
cent reduction on R. The uncertainty in the luminosity increases 
the width of possible NS radii by about 50 per cent. Another source 
of uncertainty comes from the turbulence parameter a^. With 
decreasing by an order of magnitude the spectrum hardens by 0. 1 
keV. This results in about 15 per cent decrease of the NS radius 
that is required to produce the observed spectra. Thus Qb ~ 10~ 5 
is needed to reconcile the derived NS radii with the soft equations 
of state (assuming solar composition). Such a small Qb at the same 
time yields a very large column density of the SL and a rather long 
life-time of the accreting gas in the layer (of the order of 1 s, instead 
of 10 ms as in the model of IS99). 

Finally, we would like to emphasize that our method of deter- 
mination of the NS radius from the SL spectrum is based on the 
observed color temperature of radiation alone, because the SL radi- 
ates locally at almost Eddington flux. The color temperature can be 
related to the effective temperature which is a function of the stellar 
compactness (and chemical composition) as given by equation 
This method is identical to that used for the radius-expansion X- 
ray bursts which are believed to reach Eddington luminosity (see 
e.e. lLewin. van Paradiis & Taamll993h . In contrast to the standard 
methods based on the mode ling of the thermal emission from the 
NS surface (see for example lvan Paradiis & LewiJl987t |Trumper 



2005), there is no need to know precisely either the area of the 
emitting region, or the distance to the source. 

As the standard method gives the apparent stellar radius at 
infinity, which is related to the NS parameters through 

Ro =r(i-^) 7 , (77) 

the allowed band of R and M is nearly orthogonal to that obtained 
from the color temperature and equation 1751 (see the almost ver- 
tical dashed curve in Fig. 1131 . Thus for a NS, where both the ther- 
mal emission from the surface (e.g. during the quiescence) and the 
BL emission (during the accretion phase) are observed, it would 
be possible to determine R and M independently. Interestingly our 
constraints on the N S radius are very similar to those obtained by 
Heinke et al. (2005) for the thermally emitting quiescent NS X7 in 
the globular cluster 47 Tucanae 7? = 14.5l 1 'g km. The y are also 
consis tent with the lower limit R > 14 km obtained by iTriimpeJ 
(2005) for the isolated NS RX 11856-3754. 



5 CONCLUSIONS 

We have derived the one-dimensional equations describing the SL 
model on a spherical NS surface from the usual hydrodynamic 
equations. The obtained equations are similar to those in IS99, 
except for the energy conservation law where we neglected the 
surface density of the gravitational potential energy which is of 
the second order in H/R. This difference, however, leads only 
to small quantative changes. We have also implemented a pseudo- 
Newtonian potential to account for the main general relativity cor- 
rections and considered various chemical compositions of the ac- 
creting matter. 

We have studied the vertical (radial) structure of the SL with 
different assumptions about the vertical distributions of the radi- 
ation flux and azimuthal velocity. The temperature structure and 
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the emergent radiation spectra of the SL are computed account- 
ing for the effect of Compton scattering. We showed that the lo- 
cal (at a given latitude) emergent spectra depend very little on 
details of the SL vertical structure in optically thick cases with 
Es > 100 g cm" 2 (L > 0.11/Edd)- These spectra can be described 
by the diluted Planck spectrum and are similar to the spectra of 
X-ray bursts with the same effective temperature and the effective 
surface gravity. 

The integral SL spectra were computed accounting for rela- 
tivistic effects such as the gravitational redshift and light bending, 
the relativistic Doppler effect and aberration. These spectra slightly 
depend on the inclination angle to the line of sight and on the SL 
luminosity. The local effective temperature increases with latitude, 
while the hardness factor f c decreases. This leads to only slight 
variation of the color temperature on latitude. As a result, the in- 
tegral spectra can also be well described by a single-temperature 
diluted Planck spectrum. 

We compared our theoretical integral SL spectra with the ob- 
served spectra of the LMXRBs BLs. The observed color tempera- 
tureof 2.4 ± 0.1 keV <Gilfanov et alj2003tlRevnivtsev & Gi lfanov 
2006) can be reproduced for hard equations of state of NS material. 
Our model constrains radii of NSs in LMXRBs to 13-16 km for a 
1.4 solar mass star. Soft equations of state (smaller NS radii) can 
be reconciled with the observed spectra only for very low viscosity 
Qb ~ 10 . Calculation of ab from the first principles is a chal- 
lenging problem that deserves further attention. 
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